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Abstract. The opacity of spiral galaxies is examined by 
modelling the dust and stellar content of individual galax- 
ies. The three dimensional model that we use assumes ex- 
ponential distributions for the dust and the stars in the 
disk, while the i? 1 / 4 law is used to describe the bulge. In 
this model, both absorption and scattering by the dust are 
taken into account. The model is applied to five late- type 
spiral galaxies (NGC 4013, IC 2531, UGC 1082, NGC 5529 
and NGC 5907) using their optical (and near infrared for 
IC 2531) surface photometry. For these galaxies we have 
determined the scalelengths and scaleheights of the stars 
and the dust in the disk, the bulge characteristics, the 
inclination angle and the face-on optical depth. Compu- 
tation of the dust masses, as well as the extinction as a 
function of the wavelength, are also reported. 

Having analyzed a total of seven galaxies thus far, the 
five galaxies mentioned above plus UGC 2048 and NGC 
891 presented in (Xilouris et al. 1997, 1998), we are able 
to draw some general conclusions, the most significant of 
which are: 

1) The face-on central optical depth is less than one in 
all optical bands, indicating that typical spiral galaxies 
like the ones that we have modelled would be completely 
transparent if they were to be seen face-on. 

2) The dust scaleheight is about half that of the stars, 
which means that the dust is more concentrated near the 
plane of the disk. 

3) The dust scalelength is about 1.4 times larger than that 
of the stars and the dust is more radially extended than 
the stars. 

4) The dust mass is found to be about an order of a mag- 
nitude more than previously measured using the IRAS 
fluxes, indicating the existence of a cold dust component. 
The gas-to-dust mass ratio calculated is close to the value 
derived for our Galaxy. 

5) The derived extinction law matches quite well the 



Send offprint requests to: xilouris@physics.uch.gr 



Galactic extinction law, indicating a universal dust be- 
haviour. 



Key words: galaxies: photometry - galaxies: spiral - 
galaxies: ISM - ISM: dust, extinction 

1. Introduction 

The amount and the way that dust is distributed within 
spiral galaxies as well as its extinction effects on the galac- 
tic starlight has been a matter of debate over the last few 
years. Since Disney, Davies and Phillipps (1989) proposed 
that spiral galaxies may be optically thick, a long series of 
papers has been published by different authors (see e.g. the 
introduction of Xilouris et al. 1997) with conclusions that 
often differ widely. Recent papers tend to support the idea 
of optically thin/moderately opaque spiral galaxies. This 
means that galaxies are either optically thin throughout 
their disks or that they are optically thin in the outer re- 
gions and moderately opaque at the center (e.g. Giovanelli 
et al. 1994, 1995; Kodaira & Yamashita 1996; Rdnnback 
& Shaver 1997; Domingue et al. 1998; Moriondo et al. 
1998; Kuchinski et al. 1998; Gonzalez et al. 1998). On the 
other hand, recent far infrared/submillimeter observations 
of spiral galaxies give a direct picture of the dust distri- 
bution and a good estimate of the total amount of dust 
(e.g. Alton et al. 1998a, 1998b; Haas et al. 1998; Bianchi 
et al. 1998; Krugel et al. 1998). These observations indi- 
cate that the cold dust component (~ 10 — 20 K) that is 
now detected, boosts the dust masses about an order of 
magnitude higher than previously estimated by IRAS ob- 
servations. On the other hand, a more radially extended 
dust distribution compared to that of the stars is derived. 

In our studies we perform detailed modelling of spiral 
galaxies, aiming at the determination of the parameters 
that describe the three-dimensional dust distribution in 
these galaxies. For the three dimensional distribution of 
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stars and dust we use simple and smooth geometries (e.g. 
axisymmctric exponentials for the stars and the dust in the 
disk and the i? 1 / 4 law for the bulge). There is no doubt 
that in reality spiral galaxies are clumpy with small and 
large scale structure (e.g. spiral arms with the dust trailing 
the arms in a patchy way). One has to keep in mind those 
oversimplifications. The edge-on picture though that we 
use cancels out most of the variation along the line of 
sight and allows us to use simple mathematical functions 
in order to get a mean description of the galaxies. We 
have verified this with detailed computations (Misiriotis 
et al. in preparation) which take into account the spiral 
structure. 

The two galaxies (UGC 2048 and NGC 891) that we 
have already modelled (Xilouris et al. 1997; hereafter Pa- 
per I and Xilouris et al. 1998; hereafter Paper II) us- 
ing the radiative transfer model of Kylafis & Bahcall 
(1987; hereafter KB87), show low values of the optical 
depth in the face-on orientation, indicating low opacities 
for these galaxies, while the predicted dust masses and 
the dust distributions are now supported by the far in- 
frared/submillimcter observations mentioned above. Our 
aim is to continue this kind of study with as many galaxies 
as possible in various orientations and thus improve the 
statistics. In this paper we do the analysis of the edge-on 
late-type spiral galaxies NGC 4013, IC 2531, UGC 1082, 
NGC 5529 and NGC 5907 (see Fig. 1). 

NGC 4013 is an Sbc Hubble type spiral galaxy at a 
distance of 11.6 Mpc (Bottcma 1995) and a B-band lumi- 
nosity Mb = -18.1 mag. Photometric properties of the 2D 
stellar distribution of this galaxy can be found in Bottcma 
(1995), where a disk with exponential scalelength of 2.3 
kpc and a sech 2 scaleheight of 0.7 kpc is reported to fit the 
F-band data. Barnaby et al. (1993) present bulge and disk 
decomposition of an H-band image, where an exponential 
disk with scalelength 2.5 kpc and scaleheight 0.2 kpc and a 
modified- Hubble type bulge with effective radius 0.24 kpc 
are fitted. Observations of this galaxy in the 21 cm neutral 
hydrogen line (Bottema 1995; 1996) show a highly warped 
gas disk, while the optical image of this galaxy (also pre- 
sented in Fig. 1) shows that the stellar distribution in the 
disk is not warped, at least not to the degree that should 
affect the axisymmetric modelling that we do [see also the 
optical warp curves presented in Florido et al. (1991)]. 

IC 2531 is an Sb Hubble type spiral galaxy at a dis- 
tance of 22 Mpc (Shaw et al. 1990) and M B = -18.8 mag. 
Two attempts at modelling IC 2531 have been reported 
(Wainscoat et al. 1989; Just et al. 1996). The model used 
was along the lines of KB87, but in a more simplified form, 
dealing with absorption by dust and neglecting multiple 
scattering effects. Also, only a local fit was made in sev- 
eral strips parallel to the disk's minor axis. Furthermore, 
a bulge component was not included in the model. Our 
work presents a global fit to the observed data using the 
full radiative transfer model described in KB87 and more 
realistic (bulge included) distributions. 



UGC 1082 is an Sb Hubble type spiral galaxy with 
a redshift of 0.009353 ± 0.000017 (Giovanelli & Haynes 
1993). Using the formulae given by Zombeck (1990), we 
derive a distance of 37 Mpc (assuming a value of H = 
75 km s _1 Mpc -1 ) for this object. The B luminosity is 
Mb — -18.4 mag. Visual inspection of this galaxy (see 
Fig. 1) shows an edge-on configuration of a relatively faint 
galactic disk and a bright bulge with an obvious dust lane 
passing through the disk. 

NGC 5529 is an Sc Hubble type spiral galaxy at a 
distance of 29.6 Mpc (Shaw et al. 1990) with M s = - 
19.6 mag. Visual inspection of this system shows a galaxy 
which is not quite edge-on. Despite this geometry, a dust 
lane is very well revealed and thus our axisymmetric 
model, which is able to deal with any inclination angle, 
can be applied. 

NGC 5907 is a well studied Sc galaxy at a distance of 
11 Mpc (van der Kruit & Searle 1982). It's B luminosity is 
Mb — -19.1 mag. Several photometric studies have been 
reported about this galaxy and the stellar distribution has 
been obtained in several bands. In van der Kruit & Searle 
(1982), U, J and F band data were fit to a truncated expo- 
nential disk in the R direction and a sech 2 type function 
in the z direction. The inclination angle was reported to 
be 87°, while the scalelength and scaleheight were found 
to be 5.7 kpc and 0.43 kpc respectively in the J-band. 
Barnaby and Thronson (1992) used H-band data to fit 
the stellar distribution of this galaxy. The disk was fitted 
with an exponential in the R direction with a scalelength 
of 4.0 kpc and a sech type function in the z direction with 
a scaleheight of 0.43 kpc. The bulge was well described by 
a modified Hubble profile. Morrison et al. (1994) obtained 
good quality R-band data and fitted the disk with expo- 
nential (in both directions) functions with a scalelength 
of 4.81 kpc and a scaleheight of 0.43 kpc. The inclination 
angle was reported to be 86.7°. 

In Sect. 2 we describe the observations and the data 
reduction, in Sect. 3 we briefly discuss the model and in 
Sect. 4 we present the model results for the five galaxies. A 
discussion follows in Sect. 5 with some general conclusions. 
A quick summary is presented in Sect. 6. 

2. Observations and data reduction 

Observations of the four galaxies NGC 4013, UGC 1082, 
NGC 5529 and NGC 5907 were made at Skinakas obser- 
vatory in Crete, using the 1.3 m telescope. The detector 
that was used was a Thomson 1024 x 1024 CCD camera 
with 19/im pixels, installed at the prime focus of the //7.7 
Ritchey-Cretien telescope. This arrangement corresponds 
to 0.39 arcsec/pixel and a field of view of 6.7' x 6.7'. For 
the observations of NGC 5907 we used a focal reducer 
that is installed on the telescope in order to increase the 
field of view. This arrangement produced a larger field of 
view (~ 13' x 13') and a scale of 0.75 arcsec/pixel (us- 
ing the same camera as above). The chip was run with 
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Table 1. Observational data. 



Galaxy 


Pixel size 


Observing 


seeing (") 


rms* (ma; 




Total exposure time (min) 




" (PC) 


nights 




B 


V 


I 


B 


V 


I 


NGC 4013 


0.39 (22) 


20 Jun 95 


1.3 


0.04 


0.02 


0.02 










0.39 (22) 


21 Jun 95 


1.0 


0.03 


0.02 


0.02 


140 


40 


25 


UGC 1082 


0.39 (70) 


09 Sep 94 


1.4 


0.04 


0.03 


0.03 


60 


50 


20 


NGC 5529 


0.39 (56) 


04 Aug 97 


1.2 




0.03 


0.03 




60 


40 


NGC 5907 


0.75 (40) 


07 May 97 


1.8 


0.03 


0.03 


0.02 


30 


20 


10 



* rms deviations from the least squares fit between catalogue magnitudes and calibrated magnitudes for the standard stars. 



a gain of 3.33 e~/ADU and showed a readout noise of 
5.6 e~ . Exposures were made through B, V, I broad band 
filters. The B and V passbands are comparable to those 
of Johnson's photometric system, while the I passband is 
comparable to that of Cousin's photometric system. The 
effective wavelength of the system (filter and camera) in 
each passband is 0.443^m, 0.564^m and 0.809^m for the 
B, V and I bands respectively. In order to calibrate the 
images, we used standard stars and a photometric proce- 
dure described in detail in Paper I. In Table 1 we present 
some data concerning these observations. 

The optical observations of IC 2531 were made in 
March, 1991 with the lm Australian National University 
(ANU) telescope at the Siding Spring Observatory (SSO). 
The CCD system was equipped with an EEV CCD of 
576 x 380 giving a pixel size of 0.56 arcsec at the f/8 
Cassegrain focus. Graham's (1982) standard stars were 
observed for calibration. Observations were made through 
B, V and I broad band filters with exposure times of 
1000, 600, and 500 seconds respectively. The seeing was 
estimated to be 2". The near-infrared observations of the 
same galaxy were made in February, 1995 with the 2.3 m 
ANU telescope at SSO using CASPIR (Cryogenic Array 
Spectrometer/Imager). The instrument was installed at 
the f/18 Cassegrain focus and was equipped with an SBRC 
InSb 256 x 256 detector giving a pixel size of 0.5 arcsec for 
its fast camera mode. The standard J filter and a modified 
K filter (Kn), whose long wavelength edge is tailored to 
exclude much of the thermal emission, were used for the 
imaging observations. The observation sequence included 
as many sky exposures as the galaxy in order to follow the 
fast changing sky emissions. Several standard stars were 
observed during the same night. 

3. Model 

The stellar emissivity (luminosity per unit volume) that 
we use consists of an exponential (in both radial and ver- 
tical directions) disk and a bulge described by the R 1 / 4 
law, namely 



+L b cxp(-7.67 J B 1/4 )B- 7/8 , 



(1) 



L(R, z) = L s exp 



R_ 



\z\ 



with h s and z s being the scalelength and scaleheight of 
the disk and 



B = 



^R 2 + z 2 {b/af 

R P 



(2) 



with R e being the effective radius of the bulge and a and 
b the semi-major and semi-minor axis respectively. Here 
L s and L b are the normalization constants for the stel- 
lar emissivity of the disk and the bulge respectively. The 
relations 



I s — 2L s h s , 
and 

h = 5.12L b R e , 



(3) 



(4) 



give the central value for the surface brightness of the disk 
and the bulge respectively, if the model galaxy is seen 
edge-on and there is no dust. 

For the extinction coefficient we use a double exponen- 
tial law, namely 



k\(R,z) = K A cxp 

hd z d 



(5) 



where k\ is the extinction coefficient at wavelength A at 
the center of the disk and hd and Zd are the scalelength and 
scaleheight respectively of the dust. The central optical 
depth of the model galaxy seen face-on is 



2n\z d . 



(6) 



For a more detailed description of all the parameters, 
the reader is referred to Sect. 4 of Paper I. 

The radiative transfer model that we have used is that 
described in KB87 (see also Paper I). For computational 
reasons, as in Papers I & II, the radiative transfer was per- 
formed inside a cylinder with radius R max — 3 max(/i s , h d ) 
and half height of 6 max(z s , R e ), so that practically all the 
galactic light as well as the dust are included. A Henyey- 
Greenstein phase function has been used for the scattering 
of the dust (Henyey & Greenstein 1941). The values for 
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Table 2. Global model fit parameters for NGC 4013. 



Parameter 


Units 


I band 


V band 


B band 


Is 


7 2 — 

mag/arcsec 


17.29 


±0.05 


19.08 


±0.04 


19.60 


±0.06 


Zs 


kpc 


0.20 


±0.01 


0.22 


±0.01 


0.19 


±0.01 


h a 


kpc 


1.76 


±0.1 


1.94 


±0.1 


2.55 


±0.1 


h 


mag/arcsec 


10.49 


±0.13 


11.73 


±0.11 


12.60 


±0.14 


Re 


kpc 


1.61 


±0.06 


1.42 


±0.04 


1.46 


±0.05 


b/a 




0.42 


±0.01 


0.44 


±0.01 


0.41 


±0.01 


A 




0.48 


±0.01 


0.67 


±0.01 


0.89 


±0.01 


Zd 


kpc 


0.11 


±0.01 


0.13 


±0.01 


0.13 


±0.01 


hd 


kpc 


2.60 


±0.1 


2.45 


±0.1 


2.62 


±0.3 


e 


degrees 


89.6 


±0.1 


89.7 


±0.1 


89.9 


±0.1 



Table 3. Global model fit parameters for IC 2531. 



Parameter 


Units 


K band 


J band 


I band 


V band 


B band 


/. 


mag/arcsec -3 


16.49 


±0.12 


17.45 


±0.35 


18.24 


±0.04 


19.48 


±0.05 


20.34 


±0.07 


Zs 


kpc 


0.45 


±0.02 


0.44 


±0.10 


0.43 


±0.01 


0.40 


±0.01 


0.43 


±0.01 


h s 


kpc 


5.04 


±0.1 


4.96 


±0.3 


5.05 


±0.1 


5.22 


±0.1 


6.78 


±0.1 


h 


mag/arcsec 2 


10.39 


±0.37 


10.97 


±0.40 


11.00 


±0.20 


11.85 


±0.20 


13.40 


±0.30 


Re 


kpc 


2.00 


±0.24 


1.97 


±0.10 


1.57 


±0.08 


1.23 


±0.08 


1.96 


±0.15 


b/a 




0.68 


±0.03 


0.69 


±0.02 


0.65 


±0.02 


0.65 


±0.02 


0.63 


±0.02 






0.02 


±0.01 


0.06 


±0.02 


0.22 


±0.01 


0.30 


±0.01 


0.40 


±0.01 


Zd 


kpc 


0.22 


±0.03 


0.20 


±0.07 


0.21 


±0.01 


0.23 


±0.01 


0.27 


±0.01 


h d 


kpc 


8.00 


±0.3 


8.08 


±0.6 


8.43 


±0.2 


8.18 


±0.2 


8.88 


±0.3 


9 


degrees 


89.6 


±0.2 


89.7 


±0.2 


89.6 


±0.2 


89.6 


±0.2 


89.6 


±0.2 



Table 4. Global model fit parameters for UGC 1082. 



Parameter 


Units 


I band 


V band 


B band 


I. 


mag/arcsec 2 


19.02 


±0.11 


20.50 


±0.07 


21.29 


±0.09 


Zs 


kpc 


0.44 


±0.02 


0.47 


±0.02 


0.44 


±0.03 


h 3 


kpc 


4.13 


±0.2 


4.43 


±0.1 


4.57 


±0.2 


h 


mag/arcsec 2 


9.42 


±0.08 


10.98 


±0.10 


12.17 


±0.09 


Re 


kpc 


1.34 


±0.12 


1.30 


±0.14 


1.44 


±0.16 


b/a 




0.59 


±0.01 


0.64 


±0.01 


0.63 


±0.01 






0.17 


±0.01 


0.27 


±0.02 


0.34 


±0.02 


Zd 


kpc 


0.29 


±0.02 


0.27 


±0.02 


0.27 


±0.03 


hd 


kpc 


5.47 


±0.4 


5.78 


±0.3 


5.55 


±0.4 


6 


degrees 


89.7 


±0.1 


89.7 


±0.1 


89.7 


±0.1 



the anisotropy parameter g and the albedo u> have been 
taken from Bruzual et al. (1988). Our task has been to 
find those values of the parameters in Eqs. (1) - (6) which 
create images of the model galaxies as close as possible to 
the images of the observed galaxies. For the model fitting 
techniques that we use, the reader is referred to Papers I 
and II. 



4. Results 

In Tables 2 through 6, we give the parameters of the best 
fit models to the observed data along with their 95% con- 
fidence intervals for NGC 4013, IC 2531, UGC 1082, NGC 
5529 and NGC 5907 respectively (for Table 7 describing 
NGC 891 see Sect. 5). Using the values of the parame- 
ters that describe the stellar and dust distributions of the 
galaxies, given in the above tables, we have created model 
images and have compared them with the observed images 
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Table 5. Global model fit parameters for NGC 5529. 



Parameter 


Units 


I band 


V band 


Is 


7 2 — 

mag/arcsec 


17.13 


±0.05 


18.42 


±0.04 


Zs 


kpc 


0.42 


±0.01 


0.43 


±0.01 


h s 


kpc 


4.30 


±0.1 


4.59 


±0.1 


h 


mag/arcsec 


10.98 


±0.06 


12.25 


±0.07 


Re 


kpc 


1.95 


±0.01 


1.95 


±0.01 


b/a 




0.59 


±0.01 


0.58 


±0.01 


A 




0.42 


±0.01 


0.65 


±0.01 


Zd 


kpc 


0.37 


±0.01 


0.32 


±0.01 


hd 


kpc 


7.70 


±0.2 


7.10 


±0.1 


e 


degrees 


87.2 


±0.1 


87.4 


±0.1 



Table 6. Global model fit parameters for NGC 5907. 



Parameter 


Units 


I band 


V band 


B band 


Is 


mag/arcsec 2 


17.49 


±0.02 


18.99 


±0.02 


19.86 


±0.03 


Zs 


kpc 


0.32 


±0.01 


0.34 


±0.01 


0.34 


±0.01 


h a 


kpc 


3.86 


±0.1 


4.91 


±0.1 


5.02 


±0.1 


h 


mag/arcsec 2 


9.38 


±0.04 


10.80 


±0.04 


11.88 


±0.05 


Re 


kpc 


1.02 


±0.01 


1.01 


±0.01 


0.93 


±0.01 


b/a 




0.34 


±0.01 


0.34 


±0.01 


0.45 


±0.01 






0.30 


±0.01 


0.49 


±0.01 


0.55 


±0.01 


Zd 


kpc 


0.10 


±0.01 


0.11 


±0.01 


0.13 


±0.01 


h d 


kpc 


5.30 


±0.3 


5.29 


±0.2 


5.30 


±0.2 


9 


degrees 


87.2 


±0.1 


87.2 


±0.2 


87.0 


±0.2 



Table 7. Global model fit parameters for NGC 891. 



Parameter 


Units 


K band 


J band 


I band 


V band 


B band 


/. 


mag/arcsec 2 


15.41 


±0.13 


16.01 


±0.11 


17.43 


±0.03 


18.89 


±0.05 


19.81 


±0.06 


Zs 


kpc 


0.34 


±0.01 


0.43 


±0.01 


0.38 


±0.01 


0.42 


±0.01 


0.43 


±0.01 


hs 


kpc 


3.87 


±0.1 


3.86 


±0.1 


4.93 


±0.1 


5.48 


±0.2 


5.67 


±0.2 


h 


mag/arcsec 2 


8.32 


±0.54 


9.35 


±0.16 


10.38 


±0.11 


10.97 


±0.22 


11.42 


±0.32 


Re 


kpc 


0.86 


±0.05 


0.87 


±0.05 


1.97 


±0.06 


1.51 


±0.09 


1.12 


±0.09 


b/a 




0.76 


±0.04 


0.71 


±0.08 


0.54 


±0.01 


0.54 


±0.01 


0.60 


±0.02 


H 




0.09 


±0.01 


0.21 


±0.01 


0.52 


±0.01 


0.85 


±0.01 


1.00 


±0.01 


Zd 


kpc 


0.25 


±0.03 


0.25 


±0.02 


0.24 


±0.01 


0.29 


±0.01 


0.31 


±0.02 


h d 


kpc 


8.33 


±0.4 


8.32 


±0.4 


7.54 


±0.1 


7.68 


±0.2 


8.01 


±0.3 


e 


degrees 


89.6 


±0.1 


89.6 


±0.1 


89.7 


±0.1 


89.8 


±0.2 


89.8 


±0.2 



(see Figs. 1 and 2). In Fig. 1 the I-band observations of 
all five galaxies are presented (NGC 4013, IC 2531, UGC 
1082, NGC 5529 and NGC 5907 from top to bottom). In 
Fig. 2 we show the model images of these galaxies in the 
I-band with the same scale and sequence as in Fig. 1 so 
that a direct comparison can be made between the two 
figures. In Figs. 3-7 we give a more detailed comparison 
between model and observation for each galaxy and each 
filter by showing vertical cuts along the minor axis. For 
this demonstration we use the 'folded' (around the minor 



axis) and photometrically averaged images of the galaxies, 
which are also the images that were used for the model fit 
(see Paper I). In each plot, the horizontal axis represents 
the offset (in kpc) along the vertical direction, with zero 
lying on the major axis of the disk. The vertical axis gives 
the surface brightness (in mag/arcsec 2 ). Real data are in- 
dicated by stars, while the model is shown as a solid line 
for each profile. The six profiles in each plot are vertical 
cuts at six different distances along the major axis. The 
cuts are at distances 0, 0.5/if , lh T s , 1.5/if , 2h T s and 2.5/ig and 
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Fig. 1. 1-band observations of the galaxies NGC 4013, IC 2531, 
UGC 1082, NGC 5529 and NGC 5907 (top to bottom). 

are plotted from bottom to top, with hi being the scale- 
length of the stars as derived from the I-band modelling. 
The magnitude scale corresponds to the lower profile of 
each set, the other profiles are shifted upwards (in bright- 
ness) by 2, 4, 6, 8 and 10 magnitudes respectively. Profiles 
between and 1.5/if correspond to brightness averaged 
over 4" (7.5" for NGC 5907) parallel to the major axis, 
while profiles at 2h\ and 2.5/if are averaged over 8" (15" 
for NGC 5907) along the same direction. This was done to 
reduce as much as possible the local dumpiness that may 
exist in particular areas of the galaxy. In these plots, only 
the data above the limiting sigma level (~ 3 sigma of the 
local sky) that were used for the fit are plotted, while the 
foreground stars have been removed. 

One can clearly see a good agreement between model 
and observation despite the small deviations in some 
places. It is worth remembering here that a global fit was 
done to the observed galaxy's image and the derived set 
of parameters (Tables 2-6) gives the best galaxy image 
as a whole. 

Subtracting the model image from the observation we 
obtain residual maps that give the distribution of the ab- 
solute value of the residuals throughout the galaxy's im- 
age (see Papers I & II). From these maps we can de- 
rive the percentage coverage of the galaxy's image with 
residuals below a certain value. We do so in Table 8, 
where the percentage coverage of the residuals below 
10%, 20%, 30%, 40% and 50% is given for the I-band data. 




Fig. 2. I-band models of the galaxies NGC 4013, IC 2531, 
UGC 1082, NGC 5529 and NGC 5907 (top to bottom). 



Table 8. Percentage coverage of the galaxy's image when the 
absolute value of the residuals is below a certain value. 



Galaxy 


< 10% 


< 20% 


< 30% 


< 40% 


< 50% 


NGC 4013 


46.6% 


73.2% 


87.0% 


94.3% 


97.8% 


IC 2531 


35.8% 


66.0% 


83.4% 


94.4% 


99.5% 


UGC 1082 


38.0% 


65.0% 


82.2% 


91.1% 


95.0% 


NGC 5529 


41.7% 


73.9% 


87.6% 


94.1% 


97.5% 


NGC 5907 


38.5% 


66.2% 


83.0% 


92.7% 


97.4% 



Approximately the same picture holds for the other bands 
(see Papers I & II). These numbers indicate the ability of 
our model to describe quite accurately the observed sur- 
face brightness of the galaxies that we modelled using rel- 
atively simple and smooth three dimensional distributions 
of stars and dust compared to the real ones with all the 
small scale dumpiness and spiral structure that they may 
have. 



4-1. Dust mass calculations 

Having found the parameters that best describe the dust 
distribution in each galaxy, we are able to have a good es- 
timate of the total dust mass of each galaxy (see Paper I) . 
In what follows we do these calculations and also compare 
with the total gas mass reported in the literature as well 
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Fig. 3. Vertical surface brightness profiles for NGC 4013 in the B (left), V (middle), I (right) bands. The six profiles in each 
plot are vertical cuts of the galactic plane at distances of 0, 0.5/if , lhi , 1.5/if , 2h\ and 2.5/if along the major axis as plotted 
from bottom to top respectively. Stars represent averaged profiles of the 'folded' image of the observation, while solid lines give 
the model. From bottom to top, the magnitude scale is shifted upwards (in brightness) by 0, 2, 4, 6, 8 and 10 magnitudes 
respectively. See the text for a more detailed description. 



as with the dust mass calculated using the IRAS fluxes 
(wherever measurements are available). 

For NGC 4013 the mass of the atomic hydrogen is 
reported to be 1.3 x 1O 9 M (Bottema 1995), while the 
mass of the molecular hydrogen is 8.9 x 10 8 Mq (Gomez 
& Garcia 1997), giving a total of M g = 2.2 x 1O 9 M for the 
gas mass. Using the parameters derived from our model, 
the dust mass was calculated (see Sect. 5.3 of Paper I) to 
be Md = 4.5 x 10 6 Mq. This implies a gas-to-dust mass 
ratio of 490. Using the IRAS 60^m and lOO^rn fluxes for 
this galaxy (Moshir et al. 1990) and Eq. (4) of Devereux 
& Young (1990), the dust mass is calculated to be 1.9 x 
10 6 Mq, giving a gas-to-dust mass ratio of 1200. 

In the case of IC 2531 the 21 cm line flux was measured 
to be S H , = 40.3 Jy km s" 1 (Huchtmeicr & Richter 1989) 
and thus the atomic hydrogen mass was calculated using 
Eq. (2) of Devereux & Young (1990) and is found to be 
M{Hi) = 4.6 x 10 9 Mq. We have not been able to locate 
a measurement for the 2.6 mm CO line for this galaxy in 
order to calculate the molecular hydrogen mass M(H2), 
but a crude approximation is to assume the same mass 
as that of the atomic hydrogen (see Sect. 5.3 of Paper I). 
Thus, the total gas mass is approximately M g = 9.2 x 
10 9 Mq. The dust mass for this galaxy as calculated from 
our model parameters is Md — 2.2 x 10 7 Mq, setting the 
gas-to-dust mass ratio to 420. 

For UGC 1082 the 21 cm line flux has been reported to 
be S Hl = 5.47 Jy km s" 1 (Giovanelli & Haynes 1993). This 
gives an atomic hydrogen mass of M(Hi) = 1.8 x 10 9 Mq. 
In absence of a 2.6 mm CO line measurement in order to 
calculate the molecular hydrogen mass we approximate it 
to be the same as that of the atomic hydrogen. The total 
gas mass is then M g = 3.6 x 10 9 Mq. Using the IRAS 60/xra 
and WO/im fluxes for this galaxy (Moshir et al. 1990) and 



Eq. (4) of Devereux & Young (1990), the dust mass is 
found to be 1.9 x 10 6 Mq. This implies a gas-to-dust mass 
ratio of 1900. However, using the parameters derived from 
our model, we calculate a dust mass of M d = 9.9 x 10 6 Mq 
and the gas-to-dust mass ratio then becomes 360. 

Adopting a 21 cm line flux of 40.8 Jy km s _1 (Hucht- 
meier & Richter 1989) for NGC 5529, we derive an atomic 
hydrogen mass of M(Hi) = 8.4 x 10 9 Mq. Doubling it, in 
the absence of molecular hydrogen mass information, we 
get a total gas mass of M g = 1.7 x 10 10 Mq. The 60^m 
and lOO^rn IRAS fluxes (Moshir et al. 1990) give a dust 
mass of 5.7 x 10 6 Mq which implies a gas-to-dust mass 
ratio of 3000. With our model calculations we find a dust 
mass of 3.6 x 10 7 Mq, giving a gas-to-dust mass ratio of 
470. 

For NGC 5907 exact gas mass calculations are pre- 
sented in Dumke et al. (1997). The total gas mass is re- 
ported to be M g = 7.8 x 10 9 M Q . The IRAS 60^m and 
100/iTO fluxes (Moshir et al. 1990) give a dust mass of 
3.6 x 10 6 Mq, while with our model parameters we derive 
a dust mass of 1.5 x 10 7 Mq. The gas-to-dust mass ratio 
then becomes 2200 using the IRAS measurements and 520 
using our model's parameters. 

4-2. Relative extinction values 

Adopting a mean scalchcight for the dust from all bands 
(see Sect. 5) and using the values derived for the face-on 
central optical depth (t[ — 2n\Zd) in each filter, we can 
derive the extinction coefficient n\ in each band. We have 
found that kb/kv = 1-33 and ki/kv = 0.72 for NGC 
4013, kb/kv = 1-33, ki/kv = 0.74, k,j/kv = 0.20 anad 
kk/kv — 0.06 for IC 2531, hb/hv — 1-26 and ki/kv = 
0.62 for UGC 1082, Ki /k v = 0.65 for NGC 5529 and 
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Fig. 4. Same as in Fig. 3, but for IC 2531. In addition, J and K band cuts are shown. 
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Fig. 5. Same as in Fig. 3, but for UGC 1082. 
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Fig. 6. Same as in Fig. 3, but for NGC 5529. 
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Fig. 7. Same as in Fig. 3, but for NGC 5907. 



k b /k V = 1-12 and kj / tzy = 0.61 for NGC 5907. These 
ratios are identical to the ratios of the extinction values 

Ay, I Ay. 

5. Discussion 

Having analyzed a total of seven galaxies, we are now able 
to draw some general conclusions concerning the dust and 
stellar distribution characteristics of spiral galaxies. In or- 
der to have a consistent set of parameters for all the galax- 
ies, we reanalyzed the galaxy NGC 891 (already presented 
in Paper II), but now using the i? 1 / 4 law for the bulge 
stellar distribution instead of the modified Hubble profile 
that was used in Paper II. The values derived from the 
modelling are presented in Table 7 of the present paper. 

We have quantified the opacity of the galaxies in terms 
of their central face-on optical depth . For all the galax- 
ies we have found central face- on optical depths with values 
less than one in all the bands that we analyzed, indicat- 
ing that they would appear to be transparent when seen 



face-on, despite their obvious dust lane in their almost 
edge-on orientation seen in the observations. The values 
of the optical depth in all bands are plotted in Fig. 8 for 
each galaxy. Different symbols indicate different bands. In 
particular, open circles correspond to the B-band, open 
squares to the V-band, open triangles to the I-band, solid 
squares to the J-band and solid triangles to the K-band. 

It is of interest to see how the scaleheight of the stars is 
related to that of the dust. We show this in Fig. 9, where 
the parameters z s and Zd (in kpc) as well as their ratio 
are plotted for each galaxy. The symbols have the same 
meaning as in Fig. 8. For the case of UGC 2048 (analyzed 
in Paper I), we give half of the derived scaleheight values 
since, as mentioned in Paper I, it's distance has probably 
been overestimated by a factor of 2. Doing the statistics in 
the V-band data we derive a mean ratio of z s /zd = 1.8±0.6 
implying that the dust is more concentrated to the plane 
of the disk with respect to the stars. This was of course, 
evident from the prominent dust lanes. This mean value 
of z s /zd is plotted as a horizontal dashed line. One can 
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Fig. 8. Values of the face-on central optical depth in all the 
bands that were modelled for each galaxy. Different symbols 
are for different bands. Open circles correspond to the B-band, 
open squares to the V-band, open triangles to the I-band, solid 
squares to the J-band and solid triangles to the K-band. For 
UGC 2048 the values are taken from Paper I. 



also notice that there is not a significant trend of z s with 
wavelength. This means that the main body of the disk is 
dominated by the old stellar population, at least at high z 
distances, giving roughly the same scaleheight in all bands. 

Of particular interest is the distribution of dust in the 
radial direction relative to the stars, as can be seen from 
the values that were derived for the scalelengths. In order 
to have a better view of how the scalelengths of the stars 
and the dust are related we give a plot (Fig. 10) similar 
to that of the scaleheights described above. The symbols 
here, have exactly the same meaning as that in Fig. 8. 
As with the scaleheights of UGC 2048, we adopt half the 
values of the scalelengths derived from the model. The 
statistics that were done in the V-band, give a mean ratio 
of hd/h s = 1.4 ± 0.2. This value is shown as a dashed 
horizontal line. In all seven galaxies we have found that 
the dust is radially more extended than the stars. This 
result has recently been confirmed for our Galaxy using 
COBE/DIRBE data (Davies et al. 1997a). Analysis of ISO 
data for several galaxies (Davies et al. 1997b, Alton et al. 
1998a), as well as SCUBA observations of NGC 891 in 
450/jm and 850/jm (Alton et al. 1998b) also indicate a 
more extended radial distribution for the dust than for 
the stars. Another work that also supports this, is that of 
Lcqueux & Guclin (1996) where a systematic reddening 
of background galaxies is found in the outer regions of 
the galaxy of Andromeda, suggesting the existence of dust 
even beyond the optical disk. 

It is important to notice at this point the deviations 
with wavelength shown in Figs. 9 and 10 for the dust pa- 
rameters. These deviations, which are mainly due to the 
different appearance of the clumpy nature of the galaxies 
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Fig. 9. The values of z s (top graph), Zd (middle graph) and 
z s /zd (bottom graph) for each galaxy. All lengths are given 
in kpc. For the explanation of the symbols see Fig. 8. The 
dashed line in the last graph gives the mean value of the ratio 
calculated by using the V-band data. 
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Fig. 10. The values of h s (top graph), hd (middle graph) and 
hs/hd (bottom graph) for each galaxy. All lengths are given 
in kpc. For the explanation of the symbols see Fig. 8. The 
dashed line in the last graph gives the mean value of the ratio 
calculated by using the V-band data. 
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Fig. 11. The values of M g (top graph), Md (middle graph) and 
Mg/Md (bottom graph) for each galaxy. All the masses are in 
units of Mq. For Md we give a direct comparison between 
the values as calculated by our model (filled circles) and those 
calculated using the IRAS fluxes (open circles) wherever IRAS 
data were available. The dashed line in the last graph gives the 
mean value of the gas-to-dust mass ratio. 



in different bands, give a better estimate of the real un- 
certainty which for Zd is ~ 0.05 kpc and for h<i is ~ 0.5 
kpc. 

An interesting result that comes out of our model anal- 
ysis is the dust mass calculation. For the two galaxies an- 
alyzed in the previous papers (UGC 2048; Paper I and 
NGC 891; Paper II) we find a dust mass which gives a 
gas-to-dust mass ratio very close to that derived for our 
Galaxy (the values of the dust mass for NGC 891 using the 
R 1 / 4 law profile becomes 5.6 x 10 7 Mq while the gas-to- 
dust mass ratio is now 150). About the same picture holds 
for the rest of the galaxies modelled in this paper. A bet- 
ter picture of this is given in Fig. 11, where the gas mass 
(top graph) and the dust mass (middle graph) are plotted 
(in units of Mq) for each galaxy. For the dust mass (mid- 
dle graph) we give the model calculations (filled circles) as 
well as the calculations using the IRAS data (open circles). 
The gas-to-dust mass ratio (using the dust mass derived 
from the model) is also given for each galaxy (bottom 
graph). A mean value of this ratio is M g /Md — 360 ± 160 
(close to the value of ~ 160 adopted for our Galaxy, see 
Spitzer 1978; Sodroski et al. 1994) and is plotted as a hor- 
izontal dashed line in the last graph. For all the galaxies 
we have found a dust mass about an order of a magnitude 
more than that calculated using the IRAS fluxes, which 
proves the existence of a cold dust component that has 
gone undetected by IRAS. This is evident from the middle 
graph of Fig. 11, where the two dust masses (the one cal- 
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Fig. 12. The observed (gray circles) values of A\/Ay for our 
Galaxy and the values calculated for all the galaxies analyzed. 
The symbols are explained in the plot. 



culated from our model and the other calculated using the 
IRAS fluxes) are directly compared. Our extinction model 
though, being independent of the dust temperature (both 
warm and cold dust contribute to the extinction of light), 
is able to calculate the whole amount of dust contained 
within the galaxies. Strong support to this argument (that 
vast amounts of cold dust exists in spiral galaxies) is now 
given by recent far infrared observations in longer wave- 
lengths, where the cold dust is now seen in emission (Alton 
et al. 1998a, 1998b; Bianchi et al. 1998; Haas et al. 1998; 
Kriigel et al. 1998). 

Having calculated the values of A\/Ay (see Sect. 4) 
for the galaxies modelled in this paper (NGC 891 values 
are kb/kv — 1-23, ki/kv — 0.61, kj/kv — 0.25 and 
kk/kv — 0.10 for the case where the R 1 / 4 law is used), 
we are able to show how the values of the relative ex- 
tinction in each band is compared with the values derived 
for our Galaxy (Rieke & Lebofsky 1985). This is done in 
Fig. 12, where the values of the extinction ratios A\/Ay 
are plotted as a function of the inverse wavelength (1/A). 
Each symbol here represents a different galaxy as shown 
inside the graph. One can see that the extinction law de- 
rived for all the galaxies is in very good agreement with 
that measured for our Galaxy, indicating a universal dust 
behaviour. 



6. Summary 

In this paper we have done the modelling of the dust and 
stellar content of five spiral galaxies (NGC 4013, IC 2531, 
UGC 1082, NGC 5529 and NGC 5907). With the two 
galaxies that have already been analyzed in previous work 
(UGC 2048 and NGC 891), we are able to draw some gen- 




Fig. 13. A typical model galaxy in the B-band shown at var- 
ious inclination angles (which are indicated on top of each 
model). The fainter surface brightness level at the outer part 
of the galaxy is 25 mag/arcsec 2 . 



Fig. 14. Mapping of the optical depth as a function of the 
inclination angle. White regions show the optically thick parts 
of the galaxy (with r > 1) while grey regions show the optically 
thin parts of the galaxy. 



eral conclusions and thus obtain a good picture of how a 
typical late-type spiral galaxy looks like. Assuming that 
the dust and the stellar disk follow an exponential 3D dis- 
tribution and that the bulge is described by an R 1 / 4 law 
profile, the typical parameters (derived from the mean val- 
ues) that describe the galaxy in the B-band are z s w 0.4 
kpc, h s « 5.0 kpc and I s « 20 mag/arcsec 2 for the stellar 
disk and Zd « 0.5z s and hd ~ 1.47i s for the dust disk. A 
mean B-band central face-on optical depth is « 0.8. 
For the bulge, being more dependent on the morpholog- 
ical type of the galaxy, we consider the values b/a = 0.5 
for the ellipticity, R e = 1.5 kpc for the effective radius 
and lb ~ 12 mag/arcsec 2 for the central edge-on surface 
brightness. Such a galaxy is shown in Fig. 13 in various 
inclination angles, ranging from exactly edge-on (90 °) to 
face-on (0 °). The faintest surface brightness level is 25 
mag/arcsec 2 . One thing that can clearly be seen is the ex- 
istence of a dust lane down to almost 85°. The decrease 
in the optical diameter (going from edge-on to face-on) is 
also very obvious. This is the effect of integrating along 
larger paths of light in the edge-on configuration than in 
the face-on one. In Fig. 14 we show how the "imagei" 



of the optical depth of such a galaxy varies with incli- 
nation angle. In each model-galaxy shown in this figure, 
the white regions show the parts of the galaxy that are 
optically thick (with r > 1) and grey regions show the 
optically thin parts of the galaxy. It is evident that the 
galaxy is optically thick, at least in the central regions, 
down to almost inclination angle of 60°, while the galaxy 
becomes totally transparent when viewed face-on. 
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